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Abstract 

We introduce the notion of covariance measure structure for square integrable stochas- 
tic processes. We define Wiener integral, we develop a suitable formalism for stochastic 
calculus of variations and we make Gaussian assumptions only when necessary. Our 
main examples are finite quadratric variation processes with stationary increments and 
the bifractional Brownian motion. 
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1 Introduction 



Different approaches have been used to extend the classical Ito's stochastic calculus. When 
the integrator stochastic process does not have the semimartingale property, then the pow- 
erful Ito's theory cannot be applied to integrate stochastically. Hence alternative ways have 
been then developed, essentially of two types: 

• a trajectorial approach, that mainly includes the rough paths theory (see |34] ) or the 

stochastic calculus via regularization (see |36j). 

• a Malliavin calculus (or stochastic calculus of variations) approach. 

Our main interest consists here in the second approach. Suppose that the integrator is a 
Gaussian process X = (-'^t)fg[o,T]- The Malliavin derivation can be naturally developed on 
a Gaussian space, see, e.g. [47^, [29] or [27]. A Skorohod (or divergence) integral can also 
be defined as the adjoint of the Malliavin derivative. The crucial ingredient is the canonical 
Hilbert space Ti (called also, improperly, by some authors reproducing kernel Hilbert space) 
of the Gaussian process X which is defined as the closure of the linear space generated by 
the indicator functions {l[o^t],t € [0, T]} with respect to the scalar product 

(l[0,i],l[0,s])w = ^(*,'S) (1) 

where R denotes the covariance of X. Nevertheless, this calculus remains more or less 
abstract if the structure of the elements of the Hilbert space Ti is not known. When we say 
abstract, we refer to the fact that, for example, it is difficult to characterize the processes 
which integrable with respect to X, to estimate the L^'-norms of the Skorohod integrals or 
to push further this calculus to obtain an Ito type formula. 

A particular case can be analyzed in a deeper way. We refer here to the situation 
when the covariance R can be explicitly written as 

R{t,s)= / K{t,u)K{s,u)du, 
Jo 

where K(t, s),0<s<t<Tisa deterministic kernel satisfying some regularity conditions. 
Enlarging, if need, our probability space, we can express the process X as 

Xt= f K{t,s)dWs (2) 
Jo 

where iWt)t(^[o,T] is a standard Wiener process and the above integral is understood in the 
Wiener sense. In this case, more concrete results can be proved, see [21 |9l |28]. In this 
framework the underlying Wiener process {Wt) is strongly used for developing anticipating 
calculus. In particular [28] puts emphasis on the case K{t, s) = e{t — s), when the variance 
scale of the process is as general as possible, including logarithmic scales. 
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The canonical space Ti can be written as 

7i = {K*y^ (^^([0,T])) 
where the "transfer operator" K* is defined on the set of elementary functions as 

K*{^){s) = K{T, sMs) + {ip{r) - ^(s)) K{dr, s) 

and extended (if possible) to Ti (or a set of functions contained in Ti). Consequently, 
a stochastic process u will be Skorohod integrable with respect to X if and only if K*u 
is Skorohod integrable with respect to W and J" u5X = J (K*u)5W . Depending on the 
regularity of K (in principal the Holder continuity of K and ^(t, s) are of interest) it 
becomes possible to have concrete results. 

Of course, the most studied case is the fractional Brownian motion (fBm), due to 
the multiple applications of this process in various area, like telecommunications, hydrology 
or economics. Recall that the fBm (-B/^)(g[o.T]) with Hurst parameter H € (0, 1) is defined 
as a centered Gaussian process starting from zero with covariance function 

i?(t,s) = i(t2^^ + s2^-|t-s|2^), t,SG[0,T]. (3) 

The process admits the Wiener integral representation ([2|) and the kernel K and the 
space TC can be characterized by the mean of fractional integrals and derivatives, see [21 
El Uni [Ml n\ among others. As a consequence, one can prove for any H the following Ito's 
formula 

/(Sf ) = /(O) + /* nBf)5B^ + H f nB^)s^"-^ds. 
Jo Jo 

One can also study the relation between "pathwise type" integrals and the divergence 

integral, the regularity of the Skorohod integral process or the Ito formula for indefinite 

integrals. 

As we mentioned, if the deterministic kernel K in the representation ([2|) is not ex- 
plicitly known, then the Malliavin calculus with respect to the Gaussian process X remains 
in an abstract form; and there are of course many situations when this kernel is not ex- 
plicitly known. As main example, we have in mind the case of the bifractional Brownian 
motion (hi-fBm). This process, denoted by B^'^ , is defined as a centered Gaussian process 
starting from zero with covariance 

ms) = ^{{t'^ + s^^f-\t-sr^) (4) 

where H € (0, 1) and K £ (0, 1]. When K = 1, then we have a standard fractional Brownian 
motion. 

This process was introduced in [18] and a "pathwise type" approach to stochastic 
calculus was provided in [35]. An interesting property of B^'^ consists in the expression 
of its quadratic variation (defined as usually as limit of Riemann sums, or in the sense of 
regularization, see [36]). The following properties hold true. 
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• If 2HK > 1, then the quadratic variation of B^'^ is zero. 

• If 2HK < 1 then the quadratic variation of B^'^ does not exist. 

• If 2HK = 1 then the quadratic variation of B^'^ at time t is equal to 2^~^t. 

The last property is remarkable; indeed, for HK = ^ we have a Gaussian process which 
has the same quadratic variation as the Brownian motion. Moreover, the processes is not a 
semimartingale (except for the case K = \ and H = ^), it is self-similar, has no stationary 
increments and it is a quasi-helix in the sense of J. P. Kahane [2T], that is, for all s < t, 



<2i-^|t-s|^^^. (5) 



We have no information on the form and/or the properties of the kernel of the bifr actional 
Brownian motion. As a consequence, a Malliavin calculus was not yet introduced for this 
process. On the other side, it is possible to construct a stochastic calculus of pathwise type, 
via regularization and one gets an Ito formula of the Stratonovich type (see ^35j ) 



/(i?f = /(O) + f f'iB^'^)d°B 
Jo 



H,K 

s 



for any parameters H € (0, 1) and K G (0, 1]. 

The purpose of this work is to develop a Malliavin calculus with respect to processes 
X having a covariance measure structure in sense that the covariance is the distribution 
function of a (possibly signed) measure on B{[0, T]^). This approach is particularly suitable 
for processes whose representation form ([2|) is not explicitely given. 

We will see that under this assumption, we can define suitable spaces on which the 
construction of the Malliavin derivation/Skorohod integration is coherent. 

In fact, our initial purpose is more ambitious; we start to construct a stochastic 
analysis for general (non-Gaussian) processes X having a covariance measure /i. We define 
Wiener integrals for a large enough class of deterministic functions and we define a Malliavin 
derivative and a Skorohod integral with respect to it; we can also prove certain relations 
and properties for these operators. However, if one wants to produce a consistent theory, 
then the Skorohod integral applied to deterministic integrands should coincide with the 
Wiener integral. This property is based on integration by parts on Gaussian spaces which 
is proved in Lemma 16.71 As it can be seen, that proof is completely based on the Gaussian 
character and it seems difficult to prove it for general processes. Consequently, in the 
sequel, we concentrate our study on the Gaussian case and we show various results as the 
continuity of the integral processes, the chaos expansion of local times, the relation between 
the "pathwise" and the Skorohod integrals and finally we derive the following Ito formula, 
see Corollary [8T3l for / G C'^{R) such that /" is bounded: 



f{Xt) = f{Xo) + f nx,)5X, + \ f 
Jo ^ Jo 



' nx,)d-f{s), 
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where 7(t) = Var{Xf). Our main examples include the Gaussian semimartingales, the 
fBm with H > ^, the bi-fBm with HK > ^ and processes with stationary increments. In 
the bi-fbm case, when 2HK = 1, we find a very interesting fact, that is, the bi-fBm with 
2HK = 1 satisfies the same Ito formula as the standard Wiener process, that is 

/(5f = /(O) + 1^ f\B^'^)6B^'^ + \ f"{Bf^^)ds 

where 5 denotes the Skorohod integral. 

We will also like to mention certain aspects that could be the object of a further 

study: 

• the proof of the Tanaka formula involving weighted local times; for the fBm case, this 

has been proved in [8] but the proofs necessitates the expression of the kernel K. 

• the two-parameter settings, as developed in e.g. |44| . 

• the proof of the Girsanov transform and the use of it to the study of stochastic equations 

driven by Gaussian noises, as e.g. in [3T] . 

We organized our paper as follows. In Section 2 and 3 we explain the general context 
of our study: we define the notion of covariance measure structure and we give the basic 
properties of stochastic processes with this property. Section 4 contains several examples 
of processes having covariance measure fi. Section 5 is consecrated to the construction 
of Wiener integrals for a large enough class of integrands with respect to (possibly non- 
Gaussian) process X with ^. In Section 6, for the same settings, we develop a Malliavin 
derivation and a Skorohod integration. Next, we work on a Gaussian space and our cal- 
culus assumes a more intrinsic form; we give concrete spaces of functions contained in the 
canonical Hilbert space of X and this allows us to characterize the domain of the divergence 
integral, to have Meyer inequalities and other consequences. Finally, in Section 8 we present 
the relation "pathwise" -Skorohod integrals and we derive an Ito formula; some particular 
cases are discussed in details. 

2 Preliminaries 

In this paper, a rectangle will be a subset / of of the form 

/ =]ai,6i]x]a2,&2] 

and T > will be fixed. Given F : ^ M we will denote 

A/F = F{hiM) + F{ai,a2) - F{aiM) - F{bi,a2). 

Such function will be said to vanish on the axes if F{ai,0) = F{0,a2) = for every 
ai,a2 S M+. 
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Given a continuous function F : [0,r] ^ M or a process (^t)te[o,T]! continuous in 
L^(J7), will be prolongated by convention to M by continuity. 

Definition 2.1 F : [0,T]^ — > M will he said to have a bounded planar variation if 

n 

sup \^]ta^+,]y.\t„t,+,\F\ < oo. (6) 

^ i,j=0 

where r = {0 = to < • • • < = 1} is a subdivsion of [0,T]. A function F will be said to 
be planarly increasing if for any rectangle I C [0, T]^ we have AjF > 0. 

Lemma 2.1 Let F : [0, T]^ M vanishing on the axes having a bounded planar variation. 
Then F = — F^ where F+, F^ are planarly increasing and vanishing on the axes. 

Proof: It is similar to the result of the one-parameter result, which states that a bounded 
variation function can be decomposed into the difference of two increasing functions. The 
proof of this result is written for instance in [IT] section 9-4. The proof translates into the 
planar case replacing F{b) — F{a) with AjF. ■ 

Lemma 2.2 Let F : [0, T]^ — > M_|_ be a continuous, planarly increasing function. Then 
there is a unique non-atomic, positive, finite measure n on B{[0,T]'^) such that for any 

L G ^([o,r]2) 

^(/) = AjF 

Proof: See Theorem 12.5 of [5]. ■ 

Corollary 2.3 Let F : [0,T]^ M vanishing on the axes. Suppose that F has bounded 
planar variation. Then, there is a signed, finite measure fi on I3{[0,T]'^) such that for any 
rectangle L o/ [0, T]^ 

AiF = 

Proof: It is a consequence of Lemma 12.11 and 12.21 ■ 

We recall now the notion of finite quadratic planar variation introduced in [37]. 
Definition 2.2 A function F : [0,T]^ M has finite quadratic planar variation if 

72 / {^]s^,s^+e]x]s2,S2+e]Ff dsidS2 

converges. That limit will be called the planar quadratic variation of F. 
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We introduce now some notions related to stochastic processes. Let {n,J-,P) be a 
complete probability space. Let (^t)ie[o,T] with paths in Lj^^ and (^t)te[o,T] be a cadlag 
L^-continuous process. Let t > 0. We denote by 



We set 



I^{Y,dX,t) = ^' ds 
Jo £ 

ltiY,dX,t) = fYs ^'~^'~' ds 
Jo £ 

Ce{X,Y,t) = - / {Xs+e-Xs){Ys+e-Ys)ds. 

£ Jo 

[ Yd-X (resp. [ Yd^X) 
Jo Jo 



the limit in probability of 

I-{Y,dX,t) (resp. I+{Y,dX,t)). 

Jq Yd~X (resp. Jq Yd'^X) is called (definite) forward (resp. backward) integral of Y 
with respect to X. We denote by [X, y]^ the limit in probability of Cs{X,Y,t). [X,Y]t is 
called covariation of X and Y. If X = Y, [X, X] is called quadratic variation of X, also 
denoted by [X]. 

Remark 2.4 If I is an interval with end-points a <b, then 

[ lid-X= [ lid^X = Xb- Xa. 
Jo Jo 



Let iJ^t)t&[o,T] be a filtration satisfying the usual conditions. We recall, see [36], that 
if X is an (.Ft)-semimartingale and y is a cadlag process (resp. an (.7^t)-semimartingale) 
then fl^Yd~X (resp. [y, -'^l) is the ltd integral (resp. the classical covariation). 

If X is a continuous function and y is a cadlag function then f^Yd^X coincides 
with the Lebesgue-Stieltjes integral YdX. 

3 Square integrable processes and covariance measure struc- 
ture 

In this section we will consider a cadlag zero- mean square integrable process iXt)t£[o,T] with 
covariance 

R{s,t) = Cov{Xs,Xt). 
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For simplicity we suppose that t ^ Xt is continuous in L^(J7). R defines naturally a finitely 
additive function /ir (or simply fi) on the algebra TZ of finite disjoint rectangles included in 
[0,T]'^ with values on M. We set indeed 

A typical example of square integrable processes are Gaussian processes. 

Definition 3.1 A square integrable process will be said to have a covariance measure if 

fj, extends to the Borel a-algebra ^?([0,T]^) to a signed a-finite measure. 

We recall that a{I rectangle, / C [0,r]2) = B{[0,T]'^). 

Remark 3.1 The process {Xt)te[o,T] has covariance measure if and only if R has a bounded 
planar variation, see Corollary \ 2.!^ 

Definition 3.2 Let us recall a classical notion introduced in fWj/ and 1 15^ - A process 
iXt)t£[o,T] finite energy (in the sense of discretization) if 

n-l 

sup^£;(Xi,^, -Xtj2<oo. 

Note that if X has a covariance measure then it has finite energy. Indeed for a given 
subdivision to < ti < . . . < tn, we have 

n~l n—1 n— 1 

^E{Xt^^-^ -Xt^f = ^ A]t^^t^^^]2i? < ^ \^/::^^t,,u+i\>c]t,,t,+^]R ■ 

1=0 i=0 ij=0 

Remark 3.2 Let X be a process with covariance measure. Then X has a supplementary 
property related to the energy. There is a function £ : [0, T] such that, for each 

sequence of subdivisions (r^) = {0 = to < < • • • < = T}, whose mesh converges to 
zero, the quantity 

n 

E {xti+iAt - xt^At) , (7) 

i=l 

converges uniformly in t, to £. 
Indeed 

n 

E {Xt^^.M - Xt^^t)' = KD"" n [0, tf) 

i=l 

where = Ur^o^l^i' have n?^=o ~ G [0,7"]} . From now on we will 

set 

D = {{s,s)\s G [0,T]} and Dt = D^[{),tf. 
Then, for every t E [0, T], 

£{t)=fi{Dt). 
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We will introduce the notion of energy in the sense of regularization, see [37j . 
Definition 3.3 A process (^t)£e[o,T] said to have finite energy if 

limEiCeiX,X,t)) 

e—*0 

uniformly exists. This limit will be further denoted by £{X){t). 

From now on, if we do not explicit contrary, we will essentially use the notion of 
energy in the sense of regularization. 

Lemma 3.3 If X has a covariance measure fi, then it has finite energy. Moreover 

S{X){t) = ^i{Dt). 

Proof: It holds that 

E{Ce{X,X,t)) = - f dsE {Xs+e - Xsf = - rdsA],,,+,]2fl 
^ Jo ^ Jo 

I r' r . . .1 



where 

feiyi,y2) 

We observe that 



ds dn{yi,y2) = - dfi{yi,y2)fe{yi,y2), 

£ Jo J]s,s+eP £ J[0,t+£]2 



iLe5(](yi - e) V {y2 - e) V 0, yi A ^2]) if Vi G]y2 - e, ^2 + e] 
otherwise. 



2/2)1 < - = 2 and fe{yuy2) {l [f J' 
£ t U . 1/2 7^ 2/1. 

So by Lebesgue dominated convergence theorem, 

E{Ce{X,X,t)) ^S{X){t), 
with £{X){t) = fi{Dt). m 

We recall a result established in [37j, see Proposition 3.9. 

Lemma 3.4 Let (Xt)t^[Q^T] 0, continuous, zero-mean Gaussian process with finite energy. 
Then Cs{X,X,t) converges in probability and it is deterministic for every t G [0, T] if and 
only if the planar quadratic variation of R is zero. In that case [X, X] exists and equals 
£{X). 

This allows to establish the following result. 
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Proposition 3.5 Let (^t)fg[o,T] a zero-mean continuous Gaussian process, Xq = 0, 
having a covariance measure /x. Then 

[X,X]t = fi{Dt). 

In particular the quadratic variation is deterministic. 

Proof: First, if R has bounded planar variation, then it has zero planar quadratic varia- 
tion. Indeed, by Corollary 12.31 i? has a covariance measure fi and so 



- /'dt(A]i,,+,pi?)'<- ['dt\fi\ {]t,t + e]^)T{e), 
^ Jo ^ JO 



(8) 



where 

r(e) = sup |A],^t]2i?| . 

|s-t|<e 

Since R is uniformly continuous, r(e) 0. Using the same argument as in the proof of 
Lemma [331 we conclude that ^ converges to zero. 

Second, observe that Lemma 13.31 implies that X has finite energy. Therefore the 
result follows from Lemma [37 



4 Some examples of processes with covariance measure 

4.1 X is a Gaussian martingale 

It is well known, see |37[ 140]. that [X] is deterministic. We denote A(t) = [X]t. In this case 

R{si,S2) = A(si A S2) 

so that 

fiiB)= [ 5{ds2-si)X{dsi), BeB{[0,Tf). 

JB 

If X is a classical Wiener process, then A(x) = x. The support of is the diagonal D, so ^ 
and the Lebesgue measure are mutually singular. 

4.2 X is a Gaussian (jFj-semimartingale 

We recall, see |40[ I13|. that X is a semimartingale if and only if it is a quasimartingale, i.e. 

We remark that if X is an (^f)-martingale or a process such that E {\\X\\rp) < cx), 
where H^Hj^ is the total variation, then the above condition is easily verified. According to 
[30] n extends to a measure. 
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4.3 X is a fractional Brownian motion B ,H > | 

We recall that its covariance equals, for every si,S2 G [0)^] 

R{SI,S2) = -[S^ +S2 -\S2-Si\ j. 

r\'2 r-) '2,11 2 

In that case q^^q^^ = 2H{2H — 1) |s2 — si| ~ in the sense of distributions. Since R 
vanishes on the axes, we have 

Rih,t2) = / dsi / dS2^ — TT—. 

Jo Jo dsids2 

The function R has bounded planar variation because Qg^Qg^ is non-negative. Therefore, 
for given / =]ai, 02] x]^ii &2], we have 

\AiR\ = AjR. 
Hence for a subdivision (ti)^o -^1 

TV N 

Yl \^]u,u+i]x]t,,t,+^]R\ = Yl ^]u,u+i]x]tj,t,+i]R = R{T,T). 

i,j=0 i,j=0 

So the condition ^ is verified. 



4.4 X is a bifractional Brownian motion with H G (0,1), K G (0,1] and 
2HK > 1 

We refer to [18], [35] for the definition and the basic properties of this process. The covari- 
ance of the bi-fBm is given by ^ . We can write its covariance as 

Risi,S2) = i?l(si, S2) + -^2(51,52), 

where 

Riisi, S2) = ^ [{sj^ + s^f - {sl^^ + sl^^)] (9) 

and 

R2{sr.S2) = ^ [-1.2 - s,r^ + s^^ + sr^"] . (10) 

We therefore have 

d'Rl _ 4H^K{K - 1) 2^ K-2 2H-12H-1 

Since Ri is of class C^(]0, T]^) and is always negative, -Ri is the distribution function 

of a negative absolutely continuous finite measure, having q^^qI^ for density. 

Concerning the term R2 we suppose 2HK > 1. i?2 is (up to a constant) also the 
covariance function of a fractional Brownian motion of index HK. 
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• If 2HK > 1 then R2 is the distribution ftinction of an absolutely continuous positive 
measure with density = 2HK{2HK — 1) \si — S2f^^~^ which belongs of course 

toLH[o,r]2). 

• If 2HK = 1, R2{SI, S2) = 2^(51 + S2 - |si - S2I). 



4.5 Processes with weakly stationary increments 

A process (-'^t)tG[o,T] with covariance R is said with weakly stationary increments if 
for any s,t e [0, T], h> 0, the covariance R{t + h,s + h) does not depend on h. 

Remark 4.1 // (-X't)j£[o,T] 0, Gaussian process then (-^t)te[o,T] ^^^^^ weakly stationary 

increments if and only it has stationary increments, that is, for every subdivision = to < 
ti < ... <tn and for every h>0 the law of (Xj^+ft - Xt^+h, ■■■ , ^t^+h - -'^tn-i+ft) does not 
depend on h. 

Wc consider a zero-mean continuous in process (^t)tg[o,T] such that Xq = a.s. 
Let d{s,t) be the associate canonical distance, i.e. 

d\s,t)=E{Xt-Xs)\ s,te[0,T]. 

Since {Xt)te[o,T] has stationary increments one can write 

d^{s, t) = Q{t - s), where Q{t) = d{0, t). 

Therefore the covariance function R can be expressed as 

R{s,t) = ^{Q{s)+Q{t)-Q{s-t)). 

A typical example is provided when X is a fractional Brownian motion . In that case 

Qis) = s^^. 

Remark 4.2 X has finite energy if and only if Q'{0+) exists. This follows immediately 
from the property 

EiC,{X,X,t)) = ^-^. 



We can characterize conditions on Q so that X has a covariance measure. 
Proposition 4.3 If Q" is a Radon measure, then X has a covariance measure. 

Remark 4.4 Previous assumption is equivalent to Q' being of bounded variation. In that 
case Q is absolutely continuous. 
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Proof of the Proposition 14. 3t Since 

R{S1,S2) = ^ (Q(si) + Q{S2) - Q{S2 - Si)) 

we have 

in the sense of distributions. This means in particular that for ip,ip £ P(M) (the space of 
smooth test functions with compact support) 



R{si,S2)(p'{si)'ip'{s2)dsids2 = - / ds2ipis2) / (p{si) 

JR JR 



Q{s2 - dsi] 



Example 4.5 We provide now an example of a process with stationary increments, inves- 
tigated for financial applications purposes by f^. It is called mixed fractional Brownian 
motion and it is defined as X = W + , where W is a Wiener process and is a 
fractional Brownian motion with Hurst parameter H > ^, independent from W . proves 
that X is a semimartingale if and only if H > |. 
X is a Gaussian process with 

Q(t) = |t| + (11) 

Moreover 

Q"{dt) = 26o + q'{t)dt, 

where q{t) = \t\^^~^2Hsign{t). 

The example above is still very particular. 

Suppose that Q" is a Radon measure. Then, the function Q' can be decomposed in 
the following way 

Q = Qsc + Qd + Q aci 

where Q'^^ is continuous and singular, Q^^, is absolutely continuous and Q'^ = Yln^n^^n^ 
with (x„) - sequence of nonnegative numbers and 7^ G K. 
For instance if Q is as in (llip then 

Q'sM = 0, Q'aM = 2H\t\^''-hign{t), Q'^{t) = 26o. 

A more involved example is the following. Consider Gaussian zero-mean process 
(Xj)jg[o,T] with stationary increments, Xq = such that 

Q{t) = 1 2^H-i\^^2H ; jjj - 1 
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In this case it holds that 

Q'it) 

and 



sign{t) : 1^1 ^ ^ 

2^^H\t\'^^-^sign{t) : \t\ > ^ 



Q'^ = 26o + {2H - l)6i + {2H - 1)5^ i , 

2 2 

Qac - I 2^^Hi2H - l)\t\^^-^ : \t\ > |. 

Proposition 4.6 Let (^t)fg[o,T] ^ Gaussian process with stationary increments such 
that Xq = 0. Suppose that Q" is a measure. Then 

Proof: It follows from the Proposition 13.51 and the fact that 

H{dsi,ds2) = dsiQ"{ds2 - si). 



4.6 Non-Gaussian examples 

A wide class of non-Gaussian processes having a covariance measure structure can be pro- 
vided. We will illustrate how to produce such processes living in the second Wiener chaos. 
Let us define, for every t G [0,T], 



Zt = y f{u,zi)f{u,Z2)du^ dBz^dB^^, 



where {Bz)z£R is a standard Brownian motion and / : — > M is a Borel function such that 

T f-T / r \ 2 

f{t,z)f{s,z)dz] dsdt <oo. (12) 



JO 

Now, condition (fT^ assures that belongs to L^([0,T]^). Clearly, a large class of func- 
tions / satisfies p^ . 

For example, the Rosenblatt process (see [l3]) is obtained for f{t,z) = (t — 
with A; E (i, i). In that case ^T2\i is satisfied since k > ^. 

The covariance function of the process Z is given by 

R{t,s) = j fiu,zi)f{u,Z2)dv^ (^j f{v,zi)f{v,Z2)dvJdzidz2 

/ ( / f{u,zi)f{u,Z2)f{v,zi)f{v,Z2)dzidz2]dvdu 
Jo / 
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and thus 



TT^ = / f{t,Zi)f{t,Z2)f{s,Zl)f{s,Z2)dzidZ2 

^ 2 

fit,z)f{s,z)dz 

In the case of the Rosenblatt process we get = cst.\t — 

It is also possible to construct non-continuous processes that admit a covariance 
measure structure. Let us denote by K the usual kernel of the fractional Brownian motion 
with Hurst parameter H (actually, the kernel appearing in the Wiener integral representa- 
tion ([2]) of the fBm) and consider {Nt)t£[o,T] a compensated Poisson process (see e.g. [22j). 
Then is a martingale and we can define the integral 



The covariance of Z can be written as 

i-tAs 



t 

K{t, s)dN, 







si'''] 



R{t,s)= / K{t,u)K{s,u)du = -{t'^'' + 3'^'' -\t 
Jo 2 

Then it is clear that for H > ^ the above process Z has covariance measure structure. 



5 The Wiener integral 

The Wiener integral, for integrators X which are not the classical Brownian motion, was 
considered by several authors. Among the most recent references there are [33j for the case 
of fractional Brownian motion and [20] when X is a second order process. 

We will consider in this paragraph a zero- mean, square integrable, continuous in L'^, 
process (Xt)tG[o,T] such that Xq = 0. We denote by R its covariance and we will suppose 
that X has a covariance measure denoted by ^ which is not atomic. 

We construct here a Wiener integral with respect to such a process X. Our starting 
point is the following result, see for instance [37j: if 99 is a bounded variation continuous 
real function, it is well known that 

t rt 

ipd~X = ip{t)Xt - / Xsdip,, t € [0,T]. 
Jo 

Moreover, it holds that 

lim r{s, if, dX,t) = [ ifd'X in L'^in). (13) 
'^^o Jo 
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We denote by i3y([0,T]) the space of real functions with bounded variation, defined on 
[0, T] and by C"'^([0, T]) the set of functions on [0, T] of class C^. Clearly the above relation 
(US]) holds for G C7i([0,r]). 

By S we denote the closed linear subspace of L'^{Q) generated by Jq ipd~X,ip G 
BV{[0,T]). We define $ : BVi[0,T]) ^ 5 by 

Jo 

We introduce the set as the vector space of Borel functions (^9 : [0, T] — > M such 

that 

\M'\n\-= \'P{u)\\ip{v)\d\fi\{u,v) < oo. (14) 

y[o,T]2 

We will also use the alternative notation 

ll'^'llfwi ~ / \ip <Si <p\d\n\ < oo. (15) 



[0,T] 



Remark 5.1 a) A bounded function belongs to L^, in particular if I is a real interval, 
1/ G L^. 

b) // G L^, l[o,t](/) G for any t G [0, T] . 
For if, (p £ we set 

{v>A)n= ipiu)(l){v)dfi{u,v). (16) 

Lemma 5.2 Let tp, cl) € BV{[0,T]). Then 

{if, <f>)n=E( r ifd-x r <pd-x) (17) 



and 







T 



= E( I ^d-X ] . (18) 



Proof: According to (I13p . when e — > 







E { ifis) ^'^'^ ds £ cPiu) ^"+^^ ) ■ (19) 
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converges to the right member of (|17p . We observe that p9|) equals 

\ [ dsi [ dS2^{si)<l){s2)E {{Xs,+e " {X^.+e " ^.J) 

£ JO JO 



T rT 



1 

£ JO JO 



dsi / dS2<f{si)<j){s2)l^]s^^Ss+e]x]s2,S2+e]R 



r rui i'U2 

= / dfi{ui,U2)^ / / ds2(j)is2). (20) 

J[0,T]2 e Jiui-e)+ J{u2-e)+ 

By Lebesgue dominated convergence theorem, when e ^ 0, the quantity (j20p converges to 

d/i(ui,U2)l/j(ni)</>(u2) 

[0,T]2 

and the lemma is therefore proved. ■ 
Lemma 5.3 — > (93,0)^ defines a semiscalar product on BV{[0,T]). 
Proof: The bilinearity property is obvious. On the other hand, if (/? G BV{[0,T]), 

{ip,^)^ = E(^j\d-X^ >0. (21) 

■ 

We denote by ||-||^ the associated seminorm. 

Remark 5.4 We use the terminology semiscalar product and seminorm since the property 
{(f, (/j)^ =^ 99 = does not necessarily hold. Take for instance a process 



Xt = 

where W is a classical Wiener process. 



0, t<to 
Wt-to, t>to, 



Remark 5.5 One of the difficulties in the sequel is caused by the fact that \\ ■ \\^f^^ does not 
define a norm. In particular we do not have any triangle inequality. 

Remark 5.6 If fx is a positive measure, then \\ ■ |||7-^| constitutes a true seminorm. Indeed, 
if f £ Lfj,, we have 



2 



/ \f{ui)\\f{u2)\d\l^\{uuU2) 
JfO.TP 



\f\{ui)\f\{u2)dn{ui,U2) = |||/|||L 

[o,r]2 



The triangle inequality follows easily. 
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In particular if X is a fractional Brownian motion B, H > 1/2, then ||-|||^| consti- 
tutes a norm. 

We introduce the marginal measure u associated with /x. We set 

Ki?) = |M|([o,r] xi?) 

a B eB{[0,T]). 

Lemma 5.7 If f E L^, we have 

\\f\\n<\\f\\\n\<\\f\\LHu), 

where L^{i') is the classical Hilbert space of square integrable functions on [0, T] with respect 
to u. 

Proof: The first inequality is obvious. Concerning the second one, we operate via Cauchy- 
Schwartz inequality. Indeed, 



2 

\H\ 



= / \f{ui)f{u2)\d\ii\{uuU2) 

J[0,Tf 

<\[ f\ui)d\fi\iui,U2) [ f\u2)d\ll\iuuU2)\ = r f\u)dl^{u). 
J\Q,TP J[0,T12 Jo 



'[0,T]2 V[0,T]2 

■ 

Let S be the linear subspace of constituted by the linear combinations Ojl/j, 
where li is a real interval. 

Lemma 5.8 Let u be a positive measure on B{[0,T]). Then £ (resp. C°°([0,T])) is dense 
in Lp'iv). 

Proof: 

i) Wc first prove that we can reduce to Borcl bounded functions. Let / G Lp'{u). 
We set fn = {ff\'n)\/ (— n). We have fn — > f pointwise (at each point). Consequently the 
quantity 

/ l/n - f?du 0. 
^[0,T]2 

by the dominated convergence theorem. 

ii) We can reduce to simple functions, i.e. linear combination of indicators of Borel 
sets. Indeed, any bounded Borel function / is the limit of simple functions /„, again 
pointwise. Moreover the sequence (/„) can be chosen to be bounded by |/|. 

iii) At this point we can choose f = 1b, where i? is a Borel subset of [0, T]. By 
the Radon property, for every n there is an open subset O of [0, T] with B C O, such that 
v(0\B) < This shows the existence of a sequence of /„ = lo„, where — > f in L'^{y)- 
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iv) Since every open set is a union of intervals, if / = Iqi there is a sequence of step 
functions /„ converging pointwise, monotonously to /. 

v) The problem is now reduced to take / = 1/, where / is a bounded interval. It 
is clear that / can be pointwise approximated by a sequence of C° functions /„ such that 

l/nl < 1. 

vi) Finally functions can be approximated by smooth functions via mollification; 
fn = Pn*f and (pn) is a sequence of mollifiers converging to the Dirac (5-function. 

The part concerning the density of elementary functions is contained in the previous 
proof. ■ 



We can now establish an important density proposition. 

Proposition 5.9 The set C°°([0,T]) (resp. 8) is dense in Lfj_ with respect to |HI|7^[ ^'^^ ^'^ 
particular to the seminorm 

Remark 5.10 As observed in Remark 15.51 in general ||-|||-^| does not constitute a norm. 
This is the reason, why we need to operate via Lemma \5. 7\ 

Proof of the Proposition I5.9t Let / € L^. We need to find a sequence (/„) in C°°([0,T]) 
(resp. S) so that 

ll/n ~ /lll^l ^n-^oo 0. 

The conclusion follows by Lemma 15.81 and 15.71 ■ 



Corollary 5.11 It holds that 
i) (•, •)^ is a semiscalar product on L^. 
a) The linear application 

$ : BV{[Q,T]) — > L^{Vt) 

defined by 

if — > I ifd X 
Jo 

can be continuously extended to equipped with the ||-j|^-norm. Moreover we will 
still have identity il8\) for any ip L^. 

Proof: The part i) is a direct consequence of the previous result. To check ii), it is only 
necessary to prove that $ is continuous at zero. This follows from the property (jlSp . ■ 

Definition 5.1 We will set (pdX = ^{(f) and it will be called the Wiener integral of 

ip with respect to X. 
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Remark 5.12 Consider the relation ~ on defined by 

If r-^ (t><^ - 011^ = 0. 

Denoting Lj^ as the quotient of through ~ we obtain a vector space equipped with a true 
scalar product. However is not necessarily complete and so it is not a Hilbert space. 
For the simplicity of the notation, we will still denote by its quotient with respect to the 
relation ~. Two functions (p,(p will be said to be equal in if ^ ^ 

The fact that is a metric non complete space, does not disturb the linear exten- 
sion. The important property is that L^(r2) is complete. 

Lemma 5.13 Let h be cadlag. Then 

i) J hd-X = J h^dX, 

ii) jhd+X = JhdX, 



where 



h-{u) = lim/i(s). 



Proof: We only prove point i), because the other one behaves similarly. Since h is 
bounded, we recall by Remark 15.11 that h ^ L^. We have 



^" dn= rh'dX 



with hi = - h^.du. Since 



i-T rT 

'.e u l|2 







- = / / {h^im) - h-{ui)){h''{u2) - h^{u2))d^l{ul,U2) 
Jo Jo 

and /i*" — > /i_ pointwise, the conclusion follows by Lebesgue convergence theorem. 
Corollary 5.14 If h is cadlag, then 

rT rT rT 



[ hd'X = [ hdX (= / hd+X). 
Jo Jo Jo 



Proof: Taking into account Lemma |5.13^ it is enough to show that 

rT 



[ {h- h^)dX = 0. 
Jo 



This follows because 



||/i-/i-||^= / / (h - h^){ui){h - h^){u2)dfi{ui,U2) 
J J[0,T]'2 

= ^(/i(ai) - h{ai-)){h{aj) - h{aj-))iA.{{ai,aj}) = 
and because ^ is non-atomic. 
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Remark 5.15 If I is an interval with end-points a <b, then 

J IjdX = Xk- Xa. 

This is a consequence of previous Corollary and Remark \2.4\ 



Example 5.16 The bifractional Brownian motion case: a significant subspace of 

We suppose again that HK > ^. A significant subspace included in is tlie set 
L^([0,T]) since X is a classical Brownian motion. If K = 1 and H = there is even 
equality. For K = \, and > ^, we refer to [3j. 

In the other cases, we prove that 



ll/llf^i <c(F,K,r)!|/i|,.2([o,T]). 

where C{H, K) is a constant only depending on K. It holds that 



(22) 



I •HI 



T rT 

Jo 



l/K)ll/(^2)| 
T rT 



duidu2 



{UI,U2) 



duidu2 



< C{H,K) 

+C{H,K) r r |/(ni)!|/(n2 
JO Jo 



Jo 

T rT 



f{u,)\\f{u2)\ {ul" +ul"f"\l^-'uf-'du,dU2 



\ui — U2'\^^^ '^duidu2 



A + B. 



Concerning B we refer to p]), so we have only to bound the term A. It gives 

''^ ''^ \f{m)\\f[u2)\{ul^ + ul"f-\l^-W^^-Uu2dm 

r-T f-T 



Jo 



<C{H,K) [ [ \f{ui)\\f{u2)\(,uiU2f'^''-'^+^''-'du2dni 
Jo Jo 

= C{H,K) 



'du 







<c{H,K) r \f{u)\^du r . 

Jo Jo 



L^^^-^du < C{H,K,T)\\ff 



L2([0,T])- 



Let us summarize a few points of our construction. The space given by p4p is, 
due to Remark [5TT2I a space with scalar product and it is in general incomplete. The norm 
of this space is given by the inner product (|16|) . We also define (|15|) which is not a norm in 
general but it becomes a norm when |U is a positive measure. 

We denote by H the abstract completion of L^. 
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Remark 5.17 Remark \5.1\ says that l[o,a] belongs to for any a G [0,T]. Therefore 7i 
may he seen as the closure of l\Q t]i't S [0,7"] with respect to the scalar product 

< l[0,s]> l[0,i] >H= Ris,t). 

TC is now a Hilbert space equipped with the scalar product (•, •)^; it coincides with 
()17p when restricted to L^. Ti is isomorphic to the self-reproducing kernel space. Generally 
that space is the space of v : [0,T] — > M, v{t) = E (XtZ) with Z G Therefore, if 

Z = J gdX, g G L^, we have 

^(0= / / g{s2)R{dsi,ds2). 

J[0,t] J [0,1] 

Proposition 5.18 Suppose that X is Gaussian. For f G L^, we have 

E fdX^ = E (^1^ f\s)d[X], + 2 lis,>s,)f{si)f{s2)dfl(si,S2)^ . 

Proof: It is a consequence of the Proposition 13.51 and Corollary 15.111 ■ 



6 Wiener analysis for non-Gaussian processes having a co- 
variance measure structure 

The aim of this section is to construct some framework of Malliavin calculus for stochastic 
integration related to continuous processes X, which are L^-continuous, with a covariance 
measure defined in Section 2 and Xq = 0. We denote by Co([0,T]) the set of continuous 
functions on [0, T] vanishing at zero. In this section we will also suppose that the law of 
process X on Co([0,T]) has full support, i. e. the probability that X belongs to any subset 
of Co([0, T]) is strictly positive. 

We will start with a general framework. We will define the Malliavin derivative with 
some related properties in this general, not necessarily Gaussian, framework. A Skorohod 
integral with respect to X can be defined as the adjoint of the derivative operator in some 
suitable spaces. Nevertheless, Gaussian properties are needed to go into a more practical and 
less abstract situation: for instance if one wants to exhibit concrete examples of processes 
belonging to the domain of the Skorohod integral and estimates for the norm of the 
integral. A key point, where the Gaussian nature of the process intervenes is Lemma 16.71 
We refer also to the comments following that lemma. 

We denote by Cyl the set of smooth and cylindrical random variables of the form 
F = f(^j^idX,...,jct>ndXy (23) 

where n > 1, / G C^(M") and G L^. Here J (pidX represents the Wiener integral 
introduced before Remark 15.171 
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We denote by {J^t)te[o.T] the canonical filtration associated with X fulfilling the 
usual conditions. The underlying probability space is {i},J-T,P), where P is some suitable 
probability. For our consideration, it is not restrictive to suppose that O = Cq{[0,T]), so 
that Xt{u!) = uj{t) is the canonical process. We suppose moreover that the probability 
measure P has O as support. According to Section II. 3 of [25], TC^ is dense in L^(r2), 
where 

J'C^ = {f{h,...,lm), mGN, feC^{W^), h,...,ln.en*}. 

On the other side, using similar arguments as in [23j one can prove that for every 
I G there is a sequence of random variables Zn ^ S, Zn ^ I in L'^{Q). Thus Cyl is dense 
in L'^{n). 

The derivative operator D applied to F of the form ()23p gives 



DF = Y,^^f(^ (t>idX, ...J (f>ndX 
i=l ^-^ 



In particular DF belongs a.s. to and moreover 

Recall that the classical Malliavin operator D is an unbounded linear operator from 
L^(ri) into L'^{Q;TL) where H is the abstract space defined in Section [H 

We define first the set constituted by F G L'^{Q) such that there is a sequence 

{Fn) of the form (I23p and there exists Z : Q ^ L^j^ verifying two conditions: 

i) Fn — > F in L'^{n); 

ii) E\\DF^ - Zllf^i := EJ^ {D^Fn - Z\® \D,Fn - Z\d\^i\ {u,v) ™ 0. 

The set B^'^ will be the vector subspace of L^(Q) constituted by functions F such 
that there is a sequence (F„) of the form (p3|) 

i) Fn — ^ F in L'^i^); 

ii) FpF„-DF„||^"'^°°0. 

We will denote by Z = DF the W-valued r.v. such that ^Z — DF„||^ — > 0. If 



Z £ a.s. then 



DFf^= [ Ds,FD,,Fdfi{suS2) 



'[0,T] 

Note that c B^'^ and B^'^ is a Hilbert space if equipped with the scalar 

product 

{F,G),^ = E{FG)+E{DF,DG)^. (24) 
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In general |D^'^| is not a linear space equipped with scalar product since ([T5]) is not neces- 
sarily a norm. 

Remark 6.1 Cyl is a vector algebra. Moreover, if F, G (z Cyl 

D{F ■ G) = GDF + FDG. (25) 

We prove some immediate properties of the Malliavin derivative. 
Lemma 6.2 Let F G Cyl, G G iB^'^j. Then F-G e and (2S\) still holds. 

Proof: Let (Gn) be a sequence in Cyl such that 
E{Gn — G)'^ — >„^+oo 0, 

e\ I \DGn-DG\®\DGn-DG\d\A (26) 

[J[0,T]2 J 

Since F G L°°(17), FG„ — > FG in L'^{^). Remark O implies that 

D{FGn)=GnDF + FDGn. 

So 

/ d\n\\GnDF -GDF\(E)\GnDF -GDF\ = {Cn-Cf [ \DF\\DF\d\fi\ (27) 

^[0,T]2 -'[0,T]2 

If F is of type ([23]) then 

n 

DF = Y,Zi(t>u 

i=l 

where (t)i G L^, Zi G L°°{^). Therefore the expectation of ([77]) is bounded by 



n „ 

csi. V / |(/)J•|c^|^|sf(G„-G)2ZiZJ 
.:pl■>'[o,T]2 V 



(28) 



When n converges to infinity, (|28|) converges to zero since G„ ^ G in L^(J7). On the other 
hand 

/ d\iJi\{\F{DGn-DG)\®\F{DGn-DG)\) 

J[0,T]2 

= \F\^ [ d\n\\DGn- DG\(E)\DGn- DG\. 

J[0,T]2 

Since F G L°°(r2), previous term converges to zero because of ([26|) . By additivity the result 
follows. ■ 
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We denote by L?'{Q\L^) the set of stochastic processes (?xt)jg[o,T] verifying 

E {\\u\\^\'}i\^ < oo. 

We can now define the divergence operator (or the Skorohod integral) which is an unbounded 
map defined from Dom{5) C L'^{Q.;L^) to L^(r2). We say that u € L^(J7;L^) belongs to 
Dom{6) if there is a zero-mean square integrable random variable Z € L^(r2) such that 

E{FZ) = E{{DF,u)^) (29) 

for every F ^ Cyl. In other words 



(30) 



E{FZ) =E\ Ds,Fu{s2)fJ,{dsi,ds2) ] for every F G Cyl. 

\J[0,T]^ J 

Using Riesz theorem we can see that u E Dom{5) if and only if the map 

F^E{{DF,u)^) 

is continuous linear form with respect to |Hl2,2(f7)- Since Cyl is dense in L'^{Q), Z is uniquely 
characterized. We will set ^ 

Z = I u6X. 







Z will be called the Skorohod integral of u towards X. 

Definition 6.1 //uIjq ^] € Dom{6) for any t € [0, T], we set J^UgSXs := Ugl^Q i^dXs 



Remark 6.3 // llPj] holds, then it will be valid by density for every F G 



An important preliminary result in the Malliavin calculus is the following. 

Proposition 6.4 Let u G Dom{5), F G Suppose F ■ Us6Xs G ^^(O). Then 

Fu G Dom{6) and 

cT rT 



[ F ■ Us6Xs = F I UsdXs - {DF, u) 
Jo Jo 



n ■ 



Proof: We proceed using the duality relation ([22]). Let Fq G Cyl. We need to show 



e(^Fo!^F Us6Xs - =E{{DFo 



,Fu)^). (31) 



Lemma 16.21 implies that FqF G |D ' |. The left member of ()3ip gives 
E (^FoF £ Us6Xs^ - E (Fo {DF, u)^) 

= E {{D (FoF) , u)^) - E {Fo {DF, u)^) = E{< D{FoF) - F^DF, u >„). (32) 
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This gives the right member of because of the Lemma 16.21 Remark 16.31 ahows to 
conclude. ■ 



We state a useful Fubini type theorem which allows to interchange Skorohod and 
measure theory integrals. When X is a Brownian motion and the measure theory integral 
is Lebesgue integral, then the result is stated in [29] . 



Proposition 6.5 Let {G,Q,X) be a a-finite measured space. Let u : G x [0,T] x Q 
be a measurable random field with the following properties 

i) For every x ^ G, u{x, •) G Dom{5), 
ii) 



E / dX{xi) / dX{x2) / d\fi\\u\ {xi, ■) ^ \u\ {x2, ■) < oo, 
Jg Jg ^[0,T]2 

ill) There is a measurable version inVL x G of the random field i jl^ u{x,t)6xA , 

V / x&G 

iv) It holds that 

j d\{x)E(^j^ u{x,t)6Xt^ <oo. 
Then Jq d\{x)u{x, •) G Dom{5) and 

(^j d\{x)u{x,-)^ 5Xt = j d\{x)(^j^ u{x,t)6Xt^ . 

Proof: We need to prove two properties: 
a) 

/ dX{x) \u\ (x,-) G L^{n-L^) 

JG 

b) For every F G Cyl we have 

j ^^(^)_^ u{x,t)5X^ = E (^(dF, j d\{x)u{x,-)^ ^. 

It is clear that without restriction to the generality we can suppose A to be a finite measure. 
Concerning a) we write 



E 



2 

d\{x) \u\ (x, •) 



G 



(si,S2) / dX{xi)\u\{xi,si) / d\{x2)\u\{x2,S2) 
H\) -^[o.Tja Jg jg 

(33) 

/ d\{xi)d\{x2) I d\ii\{si,S2)\u\{xi,si)\u\{x2,S2). 
JgxG J\0.T] 
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Taking the expectation of (j33|) . the result a) follows from ii). For the part b) let us consider 
F E Cyl. Classical Fubini theorem implies 



e{f j d\{x)(^j^ u{x,t)6Xt 
= j dX{x)E(^F u{x,t)5X^= j dX{x)E{{DF,u{x,-))^) 

= E < dX{x) / Ds^Fu{x,si)Ds2Fu{x,S2)dfi{si,S2)> 
[Jg ^[o,t]2 j 

= EI / dfi{si, S2)Dsj^FDs2F / dX{x)u{x,si) / dX{x)u{x,S2) 
\J\o.TP Jg Jg 



DF, / dX{x)u{x,-) 
Jg I H 

At this point the proof of the proposition is concluded. ■ 

We denote by L^,2 the set of (j) : [0, r]^ ^ R such that 
. ^{ti,-)eL^, Vti G [0,T], 

• ||0(ti,-)|||w| G Lf,. 

For (/) S -L^_2 we set 

ll'^llfw|$5|H| = / \\4>{tl,-)\\\n\U{t2,-)\\\H\d\l^\{tl,t2). 

J[0,T]2 

Similarly to [D^'^j we wih define |Di'2(L^)| and even |D1'P(L^J|, p > 2. 

We first define Cyl{L^) as the set of smooth cylindrical random elements of the 

form 

n 

Ut = Y. Mt)Gi, t G [0, T],iJe G L^, Gi G Cyl. (34) 
£=1 

On L^^2 we also define the following inner semiproduct: 



<Ui,U2>H®H= / < Ul{tl,-),U2{t2,-) >n dfl{ti,t2). 

J[0,T]2 

This inner product naturally induces a seminorm Siud we have of course 

IhlllwiK)!^! > \\u\\n«>H- 

We denote by |D"^'^(L^)| the vector space of random elements u : ^ such that 
there is a sequence (un) G Cyl{Ly) and 
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i) \\u - — >n-*oo in LP{Vt), 

ii) there \s Z : ^ — > L^^2 with H-Du^ — ^Hi-wiigiiHI — 
We convene here that 

DUn : (il,t2) > Dt^Un{t2)- 

Note that until now we did not need the Gaussian assumption on X. But this is 
essential in following result. It says that when the integrand is deterministic, the Skorohod 
integral coincide with the Wiener integral. 

Proposition 6.6 Suppose X to he a Gaussian process. Let /i G L^. Then 

T r-T 

hSXs = I hdXs. 

JQ 

Proof: Let F G Cyl. The conclusion follows from the following Lemma 16.71 and density 
arguments. ■ 



Lemma 6.7 Let F G Cyl. Then 

j-T 

Ei{DF,h)^) = E[F 



hdX^ . (35) 



Proof: We use the method given in [29], Lemma 1.1. After normalization it is possible 
to suppose that \\h\\^ = 1. There is n > 1 such that F = f hdX, J (f>idX, . . . , / 4>ndX^ , 

h, (pi, . . . , G L^, f G C^{W^). We set cpo = h and we proceed by Gram-Schmidt othogo- 
nalization. The first step is given by 



Yi = j hdX - (h,^i^^j hdX = J (pidX, 



where = /^'i^lui ^iicl so on. Therefore it is possible to find a sequence (po, . . . ,(l)n ^ L, 

h — {fi,<pi ) fill ^ 

orthonormal with respect to (•, •)^, such that 

F = f{^j <PodX, </.„dx) , / G C,~(M"+i). 

Let p be the density of the standard normal distribution in M""*"^, i.e. 
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p{x) = (2tt) 2 exp (— ^ ^127=0^1) ■ Then we have 

E{{DF,h)^) =E(^^dif(^l (PodX,...,l ct^ndXyj 



{4>i,h) 



H 



E {dof {^j <PodX, ...J (t>ndX^^ = 1^^^^ dof{y)piy)dy 
f{y)p{y)yodyo = e{f j hdx^ 



which completes the proof of the lemma. ■ 

Remark 6.8 It must be pointed out that the Gaussian property of X appears to be crucial 
in the proof of Lemma \6. 7\ Actually we used the fact that uncorrected Gaussian random 
variables are independent and also the special form of the derivative of the Gaussian kernel. 
As far as we know, there are two possible proofs of this integration by parts on the Wiener 
spaces, both using the Gaussian structure: one (that we used) presented in Nualart l29^ and 
other given in Bass |^ using a Bismut's idea and based on the Frechet form of the Malliavin 
derivative. 

7 The case of Gaussian processes with a covariance measure 
structure 

Let X = {Xt)t(^[Q^T] be a zero-mean Gaussian process such that Xq = a.s. that is contin- 
uous. A classical result of |14j (see Th. 1.3.2. and Th. 1.3.3) says that 

sup |X| G L^. (36) 

ie[o,T] 

This implies in particular that X is L^-continuous. We suppose also as in previous section 
that the law of X in Co([0, T]) has full support. 

We suppose moreover that it has covariance R with covariance measure p. Since 
X is Gaussian, according to the Section [H the canonical Hilbert space 7i ol X (called 
reproducing kernel Hilbert space by some authors) provides an abstract Wiener space and 
an abstract structure of Malliavin calculus was developed, see for instance [39^ [301 147] . 

Recently, several papers were written in relation to fractional Brownian motion and 
Volterra processes of the type Xt = /q G{t, s)dWs, where G is a deterministic kernel, see for 
instance [21 19] . In this work we remain close to the intrinsic approach based on the covariance 
as in [391 l30l |47] . However their approach is based on a version of self-reproducing kernel 
space 7i which is abstract. Our construction focuses on the linear subspace of 7i which 
is constituted by functions. 
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7.1 Properties of Malliavin derivative and divergence operator 



We introduce some elements of the Malhavin calculus with respect to X. Remark 15 . 1 71 savs 
that the abstract Hilbert space H introduced in Section 5 is the topological linear space 
generated by the indicator functions {l[o,f],t S [0,T]} with respect to the scalar product 

{^[0,t]A[0,s])'H = R{t, s). 

In general, the elements of Ti. may not be functions but distributions. This is actually the 
case of the fractional Brownian motion with H > ^, see Pipiras and Taqqu [33]. Therefore 
it is more convenient to work with one subspace of 7i that contains only functions, for 
instance L^. 

We establish here some peculiar and useful properties of Skorohod integral. 

Proposition 7.1 Let u € Cyl{L^). Then u G Dom{5) and u5X € L'p{Q) for every 
1 < p < oo. 

Proof: Let u = Gip, ip G L^, G € Cyl. Proposition 16.61 savs that ip £ Dom{6). Applying 
Proposition 16.41 with F = G and u = ip, we get that ipG belongs to Dom[5) and 

[ u5X = g[ ^5Xs- f ^{ti)Dt^GdiJL{ti,t2)- 

Jo Jo ■J[0,T]'2 

If G = g{Yi, Yn), where Yi = f cpidX, I < i < n, then 

r u6X = -y2<^j^^>'HdjgiYi,...,Yn) + g{Yi,...,Yn) C i^dX. (37) 
^0 Jo 

The right member belongs obviously to each since Yj is a Gaussian random variable and 
djQ are bounded. The final result for u € Gyl{Lfj) follows by linearity. ■ 

Remark 7.2 I^Tj ) provides an explicit expression of uSX. 



We discuss now the commutativity of the derivative and Skorohod integral. First 
we observe that if F G Cyl, (DtF) G Dom{5). Moreover, if w G Cyl{L^), {Dt^u{t2)) belongs 
to [D-'^'^(L^j^2)|- Similarly to (1-46) Ch. 1 of [29], we have the following property. 

Proposition 7.3 Let u G Cyl{L^). Then 

i-T 

/ u5X G IB^'^I 
Jo 

and we have for every t 

Dt(^J^ u6X^=ut + j^ {DtUs)5Xs. (38) 
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Proof: It is enough to write the proof for u = TpG, where G ^ Cyl oi the type 

G = g{Yi,...,Yn), 

Yi = J (j)idX, 1 < i < n. According to (j37p in the proof of Proposition 17. 11 the left member 
of ([38]) gives 

V < (t>j,i^ >H Dt{djg{Y^, ...,Yn)) + DtG C ^dX + G^(t) 
~T ^0 



Y,<'p3^^>nY. ^^5^, • • • , Yn)Ut) (39) 

j=l 1=1 

+ d,g{Yi, ...,Yn) r ijdX<p,{t) + giY,,..., YMt)- 
,=1 -^0 



3 = 

On the other hand 

n 

Dtu{s) = Hs)Y, djg{Yu . . . , y„)0, (t). 
i=i 

Applying again ([57j) . through linearity, we obtain for t € [0, T], 

r = V <Pjit) r ijdjgiYi, Yn)6X 

Jo Jo 



3 
n 



V < V >w • ■■,Yn) + djg(Yi, ...,Yn) C i^dX 

1=1 ^0 



Coming back to (|39|) we get 



i>G5X^=j^ Dt{'4jG)5X + ^{t)G. 

We can now evaluate the L^(r2) norm of the Skorohod integral. 
Proposition 7.4 Let u £ \B^'^{L^)\. Then u £ Dom{5), u5X G L^{n) and 

e( [ u6x] =E{\\u\\1^) + e( [ dii{tiM) f dfi{si,S2)Ds,ut,Dt,u, 

\Jo / \J[0,T]^ J[0,T]^ 



■ (40) 



Moreover ^ 

e(^J\6X^ <^^||n||2^l + ^^^^^d|^|(ti,t2)p.ntj|f„iy (41) 
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Remark 7.5 i) Let u,v G Polarization identity implies 

vSX^ = E{<u,v >n) 



(42) 

+E I / dfi{ti,t2) dfj.{si,S2)Ds:^utiDt2Vs2] ■ 

\y[o,T]2 ^[o,T]2 y 



ii) IfuG |D ' t/ien ul[o,f] G \B) ' {L^)\ for any t €z [0, T] and consequently to Dom{5). 

Proof (of Proposition I7.4|) : Let u G Cyl{L^). By the Proposition 17.31 since u6X G D^'^ 
we get 

f ^ uSX^ =e(^(u,D u5X^ ^ 

= £' lit^ + y Dt^UsSXs^ dfi{ti,t2)^ 

= £'(||u||^)+ / dfi{ti,t2)E iuti / Dt^UgdX, 

^ ' J\0,T]^ \ Jo 



Using again the duahty relation, we get 

rT 



+ j dfi{ti,t2) < D.ut^,Dt2U. >n^ , 



which constitutes formula (j40p . 

Moreover, using Cauchy-Schwarz, we obtain 



E(^j^ u6X^ < E (^\\uf^^ + E (^j^ ^^d\fi\{ti,t2)\\Dt,u.\\^\\D.utJ^ . (43) 

Since 



d|^l(ti,t2) ll-DutJI^ = / d\iJ.\{ti,t2)\\Dt^u\\'^ 



[0,T]2 ' ^[0,T]2 

3]) is equal or smaller than 



eI [ d\fi\{ti,t2)\\DutJl^ 



/[0,T] 

and this shows (|^T|) . 

Using the fact that Cyl{L^) is dense in |B^'^(L^)| we obtain the result. 
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7.2 Continuity of the Skorohod integral process 



It is possible to connect previous objects with the classical Wiener analysis on an abstract 
Wiener space, related to Hilbert spaces 7i, see [30| [39], 

In the classical theory the Malliavin gradient (derivative) V and the divergence 
operator 5 are well defined with its domain. For instance 5 : ©^'^(T^) — > L'^{Q) is continuous 
and D^'^('H) is contained in the classical domain. However as we said the realizations of 
u € ©^'^(W) may not be functions. 

If n G |B^'^(L^)|, it belongs to B^'^('H) and its norm is given by 



u 



1,2 



/ dn{si, S2) \\D.Us^ j|^ \\D.Us., 11^ 

J[0,T]2 J 



Classically Vu is an element of L'^{n, H^Ti), where H^Hls the Hilbert space of 
bilinear continuous functionals W (X" — > M equipped with the Hilbert-Schmidt norm. 

Given u G |B^'^(L^)| C D^'^(W), we have Du £ L'^{VL] ^^,2)- The associated gradient 
Vtt is given by 



{h,k) 



[0,T]2 



Ds-^ut^h{s2)k{t2)dfj.{si,S2)dn{ti,t2), 



where h,k € L^. Its Hilbert-Schmidt norm coincides with 



J[0,T]2 



Remark 7.6 If u £ 



UgdXs = 6{u). 



Remark 7.7 The standard Sobolev- Wiener space 3^'P(T-C), p > 2 is included in the classical 
domain of 6 and the Meyer's inequality holds: 



E\6{uW<C{p)E{\\ur^ + \\Vur^^^). 



This implies that if u £ \0^''^{L^) 



E 



/ u6X < C{p)E \ \\u\\^+ < / < Dus^,Dus2 >'H dfi{si,S2 
Jo \ l^[0,T]2 




(44) 



(45) 



Consequently this gives 



E 



[ u5X ]<C{p)e\[ \< Dus-„Dus2 >n\d\fi\{si,S2) 

Jo I l^[0,T]2 



(46) 



The last inequality can he shown in a similar way as in the case of Brownian motion. One 
applies Proposition 3.2.1 p. 158 in 129^ and then one argument in the proof of Proposition 
3.2.2 again in fS^. 
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The Meyer inequalities are very useful in order to prove the continuity of the tra- 
jectories for Skorohod integral processes. We illustrate this in the next proposition. 

Proposition 7.8 Assume that the covariance measure of the process X satisfies 

[^^{{s,t]y.{s,t\)^' <\t-s\^^^, /3>0 (47) 

for some p > 1 and consider a process u S |][!)^'^*'(L^)| such that 

[ ([ [ \DaUr\\DbUe\d\fi\{a,b)] d\fi\{r,e) <oo. (48) 
Jo \Jo Jo / 

Then the Skorohod integral process ( Zt = f/^UsdXs] admits a continuous version. 

V -"^ / te[o,T] 

Proof: We can assume that the process u is centered because the process E{us)SXs 
always admits a continuous version under our hypothesis. By ([16|). ([^7|) and ([18]) we have 

I I I I |i?an.lM(0||^fe^^^)lM(^)|%l(«,^)d|/"l(r,^)) 
JO Jo JO JO / 

<c{p)[fi{is,t]x{s,t])f~'[ [ if [ \DaUr\\DbUg\d\fi\{a,b)) d\f,\{r,9 



JO \J0 JO 



< cip,T)\t-s 



Remark 7.9 In the fBm case we have that 

Hi{s,t] X {s,t]) = 

and holds with pH > 1. In the bifractional case, it follows from J18^ that 

|^i((s,t]x(s,t])|<2i-^|t-s|2^^ 
and therefore ( [^7[ ) holds if pHK > 1. 

7.3 On local times 

We will make in this paragraph a few observations on the chaotic expansion of the local 
time of a Gaussian process X having a covariance measure structure. Our analysis is basic 
and we will only aim to anticipate a possible further study. We illustrate the fact that the 
covariance measure appears to play an important role for the existence and the regularity 
of the local time. 
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Let us use the standard way to introduce the local time L{t,x) of the process X; 
that is, for every t E [0, T] and x € M, L{t, x) is defined as the density of the occupation 
measure ^ 

HA= lA{Xs)ds, Ae 
Jo 

It can be formally written as 



L{t,x) = / 6x{Xs)ds, 
Jo 



where 5 denotes the Dirac delta function and 6x {Xs ) represents a distribution in the Watan- 
abe sense, see [37] . 

Since X is a Gaussian process, it is possible to construct related multiple Wiener-Ito 
integrals. We refer to [29| or [26] for the elements of this construction. 

There is a standard method to compute the Wiener-Ito chaos expansion of L{t,x). 
It consists in approaching the Dirac function by mean-zero Gaussian kernels Ps of variance 
e and to take the limit in L'^{Q) as e ^ 0. We get (see e.g. [12]) 



oo 



for all t G [0, T], X G M where In denotes the multiple Wiener integral of order n with 
respect to X and Hm represents the Her mite polynomial of degree m. For recent references 
on the subject one can see |30l[2l]. Note that the integral Ii is nothing else that the Wiener 
integral discussed in Section 5. One can compute the second moment of L{t,x) by using 
the isometry of multiple stochastic integrals 



m\R{s,t)^ if m = n 
lo,s]J ^™ \-[o,t]J - I if m / n 



Using standard bounds as in [_12j, it follows that the norm of ()49|) is finite if 



t rt 



.■g Jo Jo (lfi([0,ul X [0,u])||,i([0,i'] X [0,i.])|)T 

It can be seen that the existence of the local time L{t,x) as random variable in L^(Q) is 
closely related to the properties of the covariance measure /x. A possible condition to ensure 
the existence of L could be 

" r \,i[oMx[o,v])r < 

Jo {\fi{[0,u] X [0,u])\\fi{[0,v] X [0,^;])!)^ 

with P> ^. Of course, this remains rather abstract and it is interesting to be checked in 
concrete cases. We refer to [32j for the Brownian case, to [12j for the fractional Bownian 
case and to [35| for the bifractional case. 
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We also mention that the properties of the covariance measure of Gaussian processes 
are actually crucial to study sample path regularity of local times like level sets, Hausdorff 
dimension etc. in the context of the existence of local non- determinism. We refer e.g. [36] 
for a complete study of path properties of Gaussian random fields and to [4^ for the case 
of bifractional Brownian motion. 

8 Ito formula in the Gaussian case and related topics 

The next step will consist in expressing the relation between Skorohod integral and integrals 
obtained via regularization. The first result is illustrative. It does not enter into specificity 
of the examples. 

Theorem 8.1 Let (yt)te[o,T] cadlag process. We take into account the following hy- 
potheses 

a) supj<'p \ Yt\ is square integrable. 

h) Y £ |D1'2(L^)|. Moreover DY verifies 

\Dt,Yt,\<Z2, V(ii,i2) G [0,r]2 |/x|a.e. (50) 
where is a square integrable random variable. 
c) For almost all (ti,t2) G [0,T]2 

1 r^^ 

lim- / Dt.Ysds (51) 

e Jt.2-e 

exists a.s. This quantity will be denoted (Dt^Yij-)- Moreover for each s, the set oft 
such that DtYs- = DfYg is null with respect to the marginal measure v. 



c') For \fi\ almost all {ti,t2), 

1 rt2+e 

lim - / DuYsds (52) 

exists a.s. It will be denoted by {Dt^Yt2+). Moreover for each s, the set oft such that 
DfYg^ = DtYg is null with respect to measure v. 

If a), b), c) (resp. a), b), c')) are verified then Y G Dom{5) and the forward integral 
J^Yd'X (resp. the backward integral J^Yd'^X) exists and 

fYd-X= I Y6X+ [ Dt,Yt,-dpi{tiM) (53) 

JQ Jo J[0,t]2 

(resp. 

I Yd+X= I YdX+ [ Dt,Yt,+dis{ti,t2).) (54) 
Jo Jo J[0,t]^ 
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Remark 8.2 i) Condition i5U\) implies trivially 



[I \DA\[ 

\i[0,T]2 J[0, 



E\ I \Dt,Ys,\ l^^^^jDt^Ys,\j d\fi\{si,S2)d\fi\{ti,t2) <(^. (55) 
a) By Proposition \7.4\ we know that Y G Dom{5). 

Hi) Taking into account the definition if j^^YSX, it will he enough to prove h53\) and |5-^| j 
replacing t with T. 

Remark 8.3 In the case of Malliavin calculus for classical Brownian motion, see Section 
one has 

diJ,{si,S2) = S{ds2 - si)dsi. 



E{ 1 1 \DtYs\^dsdt < oo. (56) 



So f55\) becomes 

'[0,T] 



Remark 8.4 Condition f57]) (resp. (f5l|) ] may he replaced by the existence a.s. of the trace 
TrDu, where 

1 /•* 

TrDu{t) = lim - / {DY^, l].,.+,])7^ds, t G [0, T]. (57) 

This is a direct consequence of Fuhini theorem. A similar condition related to symmetric 
integral appears in JTj. 

Lemma 8.5 Let (lt)tg[o,T] ^ process fulfilling points a), h), c) of Theorem \8.1[ We set 

y/ = - f Y,ds. (58) 

e J{t-e)+ 

Then Y^ G Dom{6) and for every t 

/ Y'dX ^ / YSX in L^{n). (59) 
Jo Jo 

Proof: Again, in this proof it will be enough to set t = T. First, one can prove that if 
Y G Cyl{L^), Y' G Cyl{L^) and 

= - /' DtYrdr. (60) 

^ J s—e 

Then we can establish that Y^ is a suitable limit of elements in Cyl{L^) so that Y'^ G 
|D^'^(L^)|. We omit details at this level. Relation (|60p extends then to every Y fulfilling 
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the assumptions of the theorem. According to Proposition \7A\ € Dom{5). Relation 
(j4ip in Proposition 17.41 gives 

E ^j\Y - Y^)5X^ < E{\\Y - Y^fn) + E (^j^^^^^ dm.^t^miY,, - Y^wC^ . (61) 

We have to show that both expectations converge to zero. The first expectation 

gives 

E [ j d^,{hM){Yt,-Y^\){Yt,-Y^\)\ . (62) 

Using assumption a) of the theorem, Lebesgue dominated convergence theorem implies that 
([62]) converges to 



E 



( [ dfi{h,t2){Yt, -Yt,^){Yt, -Yt,.)] . 



For each uj a.s. the discontinuities of Y{uj) are countable. The fact that is non-atomic 
implies that previous expectation is zero. 

We discuss now the second expectation. It gives 

/ %l(ti,t2) / E{D,,{Yt,-Y,\)Ds,iYt,-Y,l))d\i,\isuS2). (63) 

J[0,T]2 -'[0,T]2 

Taking into account assumptions b), c) of the theorem, previous term converges to 



E 



if d\fi\{ti,t2) [ d\fi\{si,S2){Ds,Yt,-Ds,Yt,-){Ds,Yt,-Ds,Yt,-)] 



Using Cauchy-Schwarz this is bounded by 

E Mt) j'^ dv{s){D,Yt-D,Yt-f^ . 
This quantity is zero because of c) . ■ 

Remark 8.6 If point c') is verified (instead of c) it is possible to state a similar version of 
the lemma with Yf = ^ J^^^ Ysds. 

It is interesting to observe that convergence ([59]) holds weakly in L'^{VL) even without 
assumption c). This constitutes the following proposition. 

Proposition 8.7 Let (lt)tg[o,T] ^ process fulfilling points a), h) of Theorem \8.1[ We set 
Y^ as in 158\). Then for every t, 

/ Y^5X ^ / Y5X (64) 

weakly in L'^(Q). 
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Proof: We set again t = T. One can prove directly that belongs to Dom{6) because 
of Fubini type Proposition 16.51 Indeed, we set 

G = [0,r], v{ds) = ds, u{s,t) = Ysl]s,s+e]{i) 

and we verify the assumptions of the Proposition. Using Proposition 17.41 and points a), b) it 
is clear that -E'(/q Y^dX)"^ is bounded. Then it is possible to show that the left term in ()59p 
admits a subsequence {Jq Y^"6W) converging weakly to some square integrable random 
variable Z. 

Let F € Cyl. By duality of Skorohod integral 

^[^j^ = E {{DF,Y''")j^) 



E 



E 



( [ Ds,FY,'-f,{dsi,ds2)] 

VJ[0,T]2 J 

I / Ds,FYs,-^lidsuds2)] 



Now since X is continuous, it is not difficult to see that the 

1^1 ({si} X [0, r]) = 1^1 ([0, T] X {S2}) = 0. (65) 

Using Banach-Steinhaus theorem and the density of Cyl in L'^(Q), the convergence (j59p is 
established. For u a.s the set A^(a;) of discontinuity of Y{uj) is countable. Consequently 
|/i| ([0,T] X iV(w)) = and so 



E 



I Ds,FYs,-fiidsi,ds2)] =eI I Ds,FYs,fi{dsi,ds2)] 

= E{< DF,Y >h) = e(^ Y5X^. 



Proof of the Theorem I8.lt We set again t = T. We operate only for the forward 
integral. The backward case can be treated similarly. 
Proposition 16.41 implies that 

rT rT 

1 



r{e,Y,dX,T) =- [ dsYsf 
£ JO Jo 

= - I ds I y,l]s,s+e](i)^^t + - / dsl I dfiitut2)Dt,Ysl],^s+e]{t2) 
^ Jo Jo ^ Jo \J[0,T]^ 



h{T,e)+l2iT,s). 
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Proposition 16.51 savs that 

rT 



h{T,e)= [ 
Jo 



According to Lemma 18.51 h (T, e) converges in L^{n) to Jq Y5X. 
We observe now that l2(T, e) gives 



/ dfi{ti,t2)- [ ' dsDt.Ys. (66) 

J[0,T]2 e Jt2-e 



Assumptions b) and c) together with Lebesgue dominated convergence theorem show that 
([66l) converges in L'^{Q) to 

dfj.{ti,t2)Dtj^Yt2. 

[0,T]2 



In particular, we retrieve the result in Lemma 15.131 
Corollary 8.8 Let h be a cadlag function h : [0, T] M. Then 

T rT r-T 

hdX = hd-X = hd+X. 
Jo Jo 

Proof: This is obvious because the Malliavin derivative of h vanishes. ■. 

Corollary 8.9 Let (^t)tg[o,T] to be a process fulfilling assumptions a), b), c), c') of Theorem 
\8.1\ Then the symmetric integral ofY with respect to X is defined and 

fYd"X= frSX + l [ {Dt,Yt,+ + Dt,Yt,^)dfi{ti,t2) 
Jo Jo ^ i[0,t]2 



and 



\X,Y\t= [ {Dt,Yt,+ + Dt,Yt,^)dfi{ti,t2). 



'[0,t]2 

Example 8.10 The case of a Gaussian martingale X. 

We recall by Section [J] that [X] = A, where A is a deterministic increasing function vanishing 
at zero. Under assumption a), b), c) of Theorem 18. II 



/ Yd-X= [ Y6X+ [ d\{ti)Dt,Yt,. 
Jo Jo Jo 



Let Y be .F-progressively measurable cadlag, such that Ygd[X]s < oo a.s. In [38j it is 
also shown that f^Yd~X equals the Ito integral J^^YdX. 
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It is possible to see that the ltd integral Z = YdX verifies the duality relation 



25]) and so y G Dom{5). Moreover 



Yd-X = / Y6X. 
Jo 



We will discuss now Ito formula. Theorem 18.11 allows to state the following prelimi- 
nary formulation. 

Lemma 8.11 Let f E C^(M) such that f" is bounded. Then, for every t, f'{Xt) G Dom{6), 
and 

f{Xt) = /(Xo) + f f'{Xs)6Xs + / f"{XsMsi,S2) + \ f f"{Xs)dS{s), 

Jo J At ^ Jo 

where 

£{t)=fi{Dt), Dt = {{s,s)\s<t}, At = {{si,S2)\0<si<S2<t}. 

Proof: Ito formula for finite quadratic variation processes was established for instance by 
[37| . It says 

f{Xt) = f{Xo)+ f nX)d-X + \ f f"{X)d[X]. 
Jo ^ Jo 

Now we need to apply Theorem 18.11 For this we need to verify its hypotheses. The 
assumption a) is verified because 

sup|/(Xs)| < sup l/'l sup 

s<t t<t 

Since X is a Gaussian process, (p6]l recalls that sup^<'j. |Xj| G L^(rj). On the other 
hand, setting Yt = f'{Xt), 

Dt.Yt, = f"{Xt,)l]oMiii) 

and so b) is also verified. 

and c) is verified. Therefore 

f nX)d-X= f f'{X)5X+ [ f"{Xt,)d^{t^M). 

Jo Jo JAt 

Moreover, by Lemmas 13.31 and 13.41 we have 

^l^f"{XMX]s = lJ^f"{Xs)d£{s). 

■ 

The statement of Lemma 18.111 can be precised better. 
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Lemma 8.12 Let ^{t) = R{t,t). For every -ip € 

Jo Jo JAt 

where 

Jo ^-o7o e 

pointwise. Moreover ^ is a bounded variation function. 

Proof: The integral inside the right-hand side of (j68p gives 

* , , , R(s + e, s + e) — R(s, s) , , ^ ^ , , 
i>{s)^ ^-^ds = {h + 2l2){e,t), 







where 



/■* ds 

Ii{e,t) = / Tp{s)Aus^^]2R—, 
Jo £ 

/"* R{s + e,s) - R{s,s) 



We have 



/i(e,t)= / V(s)ii^(^s+.-^.)'-= / i^{s)d{E{Ce{X,X,s)))ds. 
Jo ^ Jo 

Taking into account Definition 13.31 and the fact that is a continuous function, 

Ii{e,t) — >e^o [ '4'{s)d£{s), 



for any t € [0,T]. It remains to control l2{s,t). 
Since 



R{s + e,s) - R{s,s) = / ho,s]x]s,s+e]dfJ'{si,S2), 

J\o,t]^ 



'[0,i]2 

we have ^ 
l2i£,t) = / dfj.{si,S2) / iPis) — , 

^[0,t]2 J]siV{s2-e)+,S2] ^ 

with the convention that ]a, 6] = if 6 < a. 
We distinguish two cases: 

• If ■S2 < si, then 

l2{e,t) = 0. 
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If S2 > Si, then 

1 



■l{j{s)ds >e^0 V'(S2)- 

]siV(s2-e),S2] 



pointwise. 
Finally 

l2{£,t) — > I dfl{si,S2)lp{s2). 



A 



and the proof of ([67|) is established. 
In particular 

-ft = £{t) + 2fi{At). 

Since £ has bounded variation and the total variation of t i — /u(At) is bounded by 
|/i|([0, T]^), then 7 is also a bounded variation function. ■ 

Corollary 8.13 Let f G C^(IR) such that f" is bounded. We set 'jt = Var{Xt). Then 

f{Xt) = f{Xo)+ f nX,)6X, + \ f r\Xs)d^{s). 
Jo ^ Jo 

Proof: It follows from Lemma 18.111 and Lemma |8 . 12 1 setting 'ip{t) = f"{Xt). ■ 

We would like to examine some particular cases. For this we decompose // into 
fi'd + A*od where for every A € i?([0, T]^) 

= fi{A n Dt), fiodiA) = h{A\Dt). 

Hence fid is concentrated on the diagonal, fiod outside the diagonal. 
We recall that 

£{t) = fi{Dt), 

where £ = £{X) is the energy function defined in Section [3j Consider the repartition 
functions Rd, Rod of fid, fJ-od- We have 

Rdisi,S2) = fJ-d (]0, si] x]0, S2]) = £{si A S2) 

and 

Rod{si, S2) = flod (]0, Si] x]0, S2]) . 

Remark 8.14 i) Setting ip = £, there is a Gaussian martingale M with covariance Rd- It 
is enough to take Mt = where {Wt) is a classical Brownian motion, 

ii) Rod{si,S2) = Cov{Xs„Xs^) - Cov{Ms„Ms,). 
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Proposition 8.15 Suppose that 

i) 8 is absolutely continuous, i.e. there is a locally integrable function £' such that 

£{t) = [ £'{s)ds : 

ii) fiod is absolutely continuous with respect to Lebesgue measure. In particular one 

Aiod(]0,si]x]0,S2]) = / dsids2. 

J]0,si]x]0,S2] 0SlOS2 

In fact it is clear that 

o = r, on 0,T \Dt. 
OSIOS2 OSIOS2 

Then the conclusion of the Theorem \8.1\ holds replacing assumptions 

• Y is cadlag, 

• c) (resp. c'), 
with 

• For t a.e. Lebesgue 

1 /■* 

DtYt- = lim — / DfYgds exists a.s. 

t-s 



e^O £ 

(resp. 

1 



i+e 



DfYt-f- = lim — / DfYgds exists a.s.) 
Moreover the conclusion of the theorem can be stated as 

/ Yd-X= Y6X+ DtYt-£'{t)dt+ Dt,Yt,—^{ti,t2)dti,dt2. 
Jo Jo Jo J\0,T]2 dsids2 



[0,T]2 

(resp. 

cT rT rT r q^r 



[ Yd+X= [ Y6X+ [ DtYt+£'{t)dt+ [ Dt.Yt,^^ 
Jo Jo Jo J[0,T]^ dsids2 



{ti,t2)dti,dt2.) 
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Remark 8.16 If c) and c') of Theorem \8. 11 with Il6y\) and ( |7^[ j are verified then 
and 



[ Yd°X= [ Y6X+ [ {D,Ys+ + D,Ys^)£'{s)dt+ [ Dt.Yt,^ 
Jo Jo Jo J[o,ti2 (ysi 



ds2 



{ti,t2)dtidt2 



[X, Y]t = [\d,Ys+ + D,Ys-)£'{s)dt. 
Jo 



Before proceeding to the proof, we recall a basic result which can be found in 
Lemma 8.17 Let g G 1 < p < oo. We set 

1 rx+e 

9e{x) = - g{y)dy or ge{x) = - g{y)dy. 

^ J x—e ^ J X 

Then ge ^ g a.e. and in L^. 



Proof of Proposition I8.15t The proof follows the same line as the proof of Theorem 18.11 
a) First we need to adapt Lemma [531 to show that lim^^Q Y^dXs = J^Ys6Xs in L^(i7) 
where Y^ still denotes the same approximation process. Again we need to show that the 
right hand side of (j6T]) converges to zero when e — > 0. Its first term gives (/i + /2)(e), where 



Lemma 18.171 implies that Y^ — > Y a.e. dP ® Leb. Lebesgue dominated convergence 
theorem and Assumption a) imply that /i(e) — > and /2(£) — ^ 0, when e converges to 
zero. It remains to control the second term in (|6ip which is given by (|63p . This second 
term gives Ki{e) + K2{e) with 



dsids2 

T 



Ki{e) = E i / dsids2 

yj[o,T]2 

K2{e) = E ds\£'{s)\ £ dtiDsYf - DsYtf)^ 



\D,,Y,[-Ds,Yt,\\D,,Yl-Ds,Yt,\ 



Point b) and ([M]) allow to show that Ki{e) + K2{e) — > 0. 

b) The other point concerns the convergence of l2{T,e) appearing in the proof of 
Theorem 8.1. To prove the convergence of (|66|) we separate again ^ = fi^ + and we 
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use on the diagonal. Finally Lemma [5TT71 ([^Uj) and Lebesgue dominated convergence 
theorem show that for ti, t2, ti ^ t2 a.e. 



1 n 

' dsAi^^s ^ Alicia a.e. dP®dtidt2. 

t2 



e 



Example 8.18 Let us apply the obtained results to some particular examples. 

a) Case of a Gaussian martingale with absolutely continuous quadratic variation X{t) 
X{0) + JiX{s)ds. 

R{tut2) = A(ti A t2), 7(*) = Var{Xt) = A(t), 
S'{t) = X{t), R = Rd, 

— 7: — = a.e. Lebesgue. 
OT1OT2 

Let Y be as in Proposition 18.151 Then 

f-T f-T rT 



[ Yd-X= [ Y6X + [ DtYt-X{t)dt. 
Jo Jo Jo 



b) The case of fractional Brownian motion H > 1/2. 
We have 

R = Rod 

d^R 



^^^^^^ 2H{2H-l)\t2-t^\'^-\ i{t)=t'^. (73) 
One obtains the classical results for fractional Brownian motion as in [3], for instance 



/ Yd-X= [ Y5X + H{2H -I) [ Dt^Yt^\t2 - ti 

Jo Jo J[0,T]^ 



'^"-^dtidt2. 



'[0,T]2 

Corollarv 18.131 provides the following Ito formula: 



f{Xt) = /(Xo) + f f{X)6X + H f f"{X,)s'''-'ds. 
Jo Jo 



(74) 



c) The case of bifractional Brownian motion X = B^'^ , HK > ^. 
It is easy to verify that Var{Xt) = t^^^ so that 



Corollary 18.131 implies 

f{Xt) = f{0)+ f f'{Xs)SXs + HK f f"{Xs)s^''''-'ds. 
Jo Jo 
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In particular, if HK = 2 we get a formula which looks very similar to the one for the 
Brownian motion: 

f{Xt) = f{Xo) + f{Xs)5Xs + f"{Xs)ds. 
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